DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 93943-9008 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 



STUDY OF OFF-AXIS RADIATION EXPOSURE FROM 
RELATIVISTIC ELECTRONS TRAVERSING 
THROUGH MATTER 

by 

Richard D. Fitzpatrick 
December 1985 

Thesis Advisor: X.K. Maruyama 

Approved for public release: distribution is unlimited. 

T226304 



UNCLASSIFIED 



iURITY CLASSIFICATION OP THIS PAGE 



REPORT DOCUMENTATION PAGE 



report security classification 



1b. RESTRICTIVE MARKINGS 



SECURITY CLASSIFICATION AUTHORITY 



DECLASSIFICATION / DOWNGRADING SCHEDULE 



3 DISTRIBUTION /AVAILABILITY OF REPORT 

Approved for public release 5 
distribution is unlimited. 



PERFORMING ORGANIZATION REPORT NUMBER(S) 



S. MONITORING ORGANIZATION REPORT NUMBER(S) 



NAME OF PERFORMING ORGANIZATION 

aval Postgraduate School 



6 b OFFICE SYMBOL 
(If applicable) 

Code 61 



7a. NAME OF MONITORING ORGANIZATION 

Naval Postgraduate School 



ADDRESS (City, State , and ZIP Code) 

onterey, California 93943-5100 



7b. ADDRESS (City, State, and ZIP Code) 

Monterey, California 93943-5100 



NAME OF FUNDING/SPONSORING 
ORGANIZATION 



8 b. OFFICE SYMBOL 
(If applicable) 



9. PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 



ADDRESS (City, State, and ZIP Code) 



10. SOURCE OF FUNDING NUMBERS 



PROGRAM 


PROJECT 




TASK 


WORK UNIT 5 


ELEMENT NO 


NO 




NO 


ACCESSION NO | 
1 



TITLE (Include Security Classification) 

UDY OF OFF-AXIS RADIATION EXPOSURE FROM RELATIVISTIC ELECTRONS 
AVERS ING THROUGH MATTER 



PERSONAL AUTHOR(S) 

itzpatrick, Richard D. 



1 TYPE OF REPORT 

aster's Thesis 


13b TIME COVERED 
FROM TO 


14 DATE OF REPORT (Year, Month, Day) 

1985, December 


15 PAGE COUNT | 

144 


SUPPLEMENTARY NOTATION 



COSATI CODES 


FIELD 


GROUP 


SUB-GROUP 















18 SUBJECT TERMS {Continue on reverse if necessary and identify by block number) 

Off-Axis Radiation; Dosimetry; CYLTRAN 



BSTRACT ( Continue on reverse if necessary and identify by block number) 

Results are presented from the study of the off-axis radiation field 
caused by a relativistic, mono-directional electron beam passing through 
water . Using the Naval Postgraduate School 100 MeV linear accelerator, 
°fi~sxis radiation dose was measured with calcium fluoride thermo- 
luminescent dosimeters placed at various path lengths out to two 
radiation lengths. The off-axis dose was calculated using the electron 
transport code CYLTRAN of the integrated TIGER series of coupled 
elect ron/photon Monte Carlo transport codes .Calculations were performed 
at Los Alamos National Laboratory. Comparison with the results is 
made and CYLTRAN is found to be in agreement with experimentally 
measured values. The extension of results from one med ium ( water ) to 
another (air) appears to be valid. 



J 



DISTRIBUTION / AVAILABILITY OF ABSTRACT 
3 UNCLASSIFIED/UNLIMITED □ SAME AS RPT □ DTIC USERS 



21 



ABSTRACT SECURITY CLASSIFICATION 

Unclassified 



NAME OF RESPONSIBLE INDIVIDUAL 

Prof. Fred R. Buskirk 



22b TELEPHONE (Include Area Code) 

(408) 646-2765 



22c OFFICE SYMBOL 

Code 61Bs 



) FORM 1473, 34 MAR 



83 APR edition may be used until exhausted 
All other editions are obsolete 
1 



SECURITY CLASSIFICATION OF THIS p AGE 



.] 



Approved for public release; distribution is unlimited. 



Study of Off-Axis Radiation Exposure From Relativistic Electrons 
Traversing Through Matter 

by 

Richard Douglas Fitzpatrick 
Lieutenant, Unite# States Navy 
B.S., United States Naval Academy, 1979 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN PHYSICS 



from the 



NAVAL POSTGRADUATE SCHOOL 
December 1985 



ABSTRACT 



Results are presented from the study of the off- axis radiation field 
caused by a relativistic, mono-directional electron beam passing through 
water. Using the Naval Postgraduate School 100 MeV linear accelerator, 
off-axis radiation dose was measured with calcium fluoride thermo- 
luminescent dosimeters placed at various path lengths out to two radiation 
lengths. The off-axis dose was calculated using the electron transport code 
CYLTRAN of the integrated TIGER series of coupled electron/photon Monte 
Carlo transport codes. Calculations were performed at Los Alamos National 
Laboratory. Comparison with the results is made and CYLTRAN is found to be 
in agreement with experimentally measured values. The extension of 
experimental results in one medium (water) to another (air) appears to be 

valid. 



3 



TABLE OF CONTENTS 



I. INTRODUCTION 6 

II. BACKGROUND 8 

III. ISSUES 13 

A. EXPERIMENTAL GEOMETRY 13 

B. ENERGY DEPENDENCE M 

C. COMPUTATION M 

D. DOSIMETRY 15 

IV. CALCULATIONS 17 

V. EXPERIMENTS 31 

VI. CONCLUSIONS 48 

LIST OF REFERENCES 69 

APPENDIX A: DESCRIPTION OF CYLTRAN AND OUTPUT 70 

APPENDIX B: CYLTRAN OUTPUT OF 1 00 MeV ELECTRONS IN WATER 97 

APPENDIX C: CYLTRAN OUTPUT OF 100 MeV ELECTRONS IN AIR 107 

APPENDIX D: CYLTRAN OUTPUT OF 100 MeV ELECTRONS IN LN 2 116 

APPENDIX E: CYLTRAN OUTPUT OF 20 MeV ELECTRONS IN WATER 126 

APPENDIX F: EQUATIONS USED FOR CURVE FITTING 135 

APPENDIX G: EXPERIMENTAL RESULTS OF 100 MeV ELECTRONS IN WATER 139 

INITIAL DISTRIBUTION LIST 142 



4 



ACKNOWLEDGMENT 



Special thanks are due to Dr. Joseph Mack of Los Alamos National 
Laboratory and Ms. Louise Miles of the Naval Surface Weapons Center who 
both contibuted significant time and expertise towards this project. Also 
deserving recognition is Mr. Don Snyder of the Naval Postgraduate School 
Accelerator Laboratory whose involvement with the students proved 
invaluable. 



5 



I. INTRODUCTION 



Recent development of new high energy, high current accelerators has 
renewed interest in the ability to accurately predict radiation exposure in the 
vicinity of a monodirectional electron beam. In addition to personnel safety 
requirements, the areas of vulnerability and lethality as they apply to 
charged particle beams are becoming increasingly important. The majority 
of the literature on radiation dosimetry applies to low or moderate energy 
sources. The experimental verification of the calculations as they apply to 
high energy electron beams has been sparse and comparison of several 
computational methods has been inconclusive. 

Recent experiments at the Naval Postgraduate School by P.F. Cromar 
attempted to address some of these questions by studying the off-axis 
radiation field in liquid nitrogen as a function of radiation length [Ref. 1: p.4]. 

Focused electrons of 100 MeV incident energy were delivered to insulated 
containers filled with liquid nitrogen. The measured dose was compared to 
calculations using the electron transport code CYLTRAN . Cromar's results 
compared favorably (within a factor of 2) with the results predicted by the 
calculations. However, where comparisons could be made, Cromar's 

results appeared to differ from an earlier calculation by R.A. Lindgren 
using the code ETRAN-16, which referred to 50 Mev electrons in air 
[Ref. I: p.39] . It is the objective of this thesis to address, through a series 
of experiments, several questions concerning the validity of calculations and 
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experimental methods as they apply to high energy electron beams. Passive 
dosimetry will be used to determine: 



1) Can the computer computations accurately predict the experimental 
results? 

2) Can the results in one medium be used to predict the response in a 
second medium? 

3) Is there an energy dependence to the response? 

4) If the responses are different, what is the significance as it applies to 
the energy spectrum? 
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II. BACKGROUND 



A. INTERACTION OF RADIATION WITH MATTER 

When a primary electron passes through an absorber it will only have a 
well-defined range for energies of less than a few MeV. At these low 
energies, energy loss is by ionization and excitation of the atomic electrons 
in the atoms of the absorber[Ref.2: p . 508] . The energy loss of electrons due 
to ionization is given by [Ref. 3: p. 38]: 



I'dE'j _ 4tf8 4 n 
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where e = charge on the electron (cgs) 
n = electron density 
m = mass of electron 
c = speed of light 
I = ionization potential 
13 = ratio of electron velocity to c 



As the electron energy increases, radiative energy loss becomes more 
important than ionization. This energy loss is given by [Ref. 3: p . 68]: 
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where r 0 = classical radius of the electron 

h = Planck's constant 
Z = atomic number 

u = frequency of the highest energy 
photon produced 

The energy loss by radiation, called bremsstrahlung, is proportional to 
Z 2 of the material and increases linearly with the electrons energy, whereas 
losses due to ionization and excitation are only proportional to Z. Thus the 
radiation loss predominates at the higher energies. The approximate ratio 
of the two losses is [Ref. 3: p . 64] 
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There is a critical energy, E c , for which the two energy losses are equal: 



If, in a given material, the incident energy is much larger than the 
critical energy, then radiation loss predominates over ionization loss and can 
be considered to be the only significant loss mechanism [Ref. 3: p . 72 ] . 
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When electrons of high initial energy pass through matter, they lose 
their energy mainly by emitting bremsstrahlung photons. These high-energy 
photons, or gamma rays, in turn lose their energy mainly by the electron 
pair-production process. By a series of similar interactions and a 
succession of these simple elementary processes, a multiplication process 
comes about. An avalanche of electrons and photons develops and is referred 
to as cascade shower or electron-photon cascade. [Ref .4: p. 91 ][Ref . 5: p. 22 1 ] 
The energy of the primary electron is divided among the cascade 
particles, their sum equaling the primary energy except for the amount lost 
due to ionization. A rapid multiplication of shower particles results with the 
number of particles growing exponetially as long as their energies exceed the 

critical energy, E c . As the number of particles reaches a maximum value 

and the energy per particle is smaller, the ionization losses become more 
important than the bremsstrahlung losses. Multiplication of shower 
particles ceases and as energy is spent in ionization and excitation processes 
the number of particles decreases. [Ref. 4:pp. 92-93] 

The maximum number of shower particles amounts to : 

_£o 

^ max 

This maximum is observed at an absorber depth of: 
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t In ( Eo/ Ec ) 

L max - 

In 2 

with t max being measured in radiation lengths. [Ref. 4: p.93] 

Photon production by bremsstrahlung and their angular distribution is 
largely a statistical process. The spectrum of energies produced by an 
electron beam will be broad and the percentage of energy transferred to any 
given photon following an interaction will follow a probability distribution 
function determined by the energy of the incident electron and the atomic 
number of the material. The result is that many secondary electrons which 
are emitted cause further ionization and excitation, and may penetrate far 
beyond the range of the primary electrons [Ref. 2: p . 509] . The precise 

mathematical development of showers is complicated and can contain large 
fluctuations from the average behavior. In practice, it becomes necessary 
to use Monte Carlo techniques to calculate the shower induced dose. 

The difficulty in separating the effects of the primary and secondary 
electrons makes range an inappropriate measure of electron absorption 

above a few MeV energy [Ref. 2: p . 509 . ] . The concept of radiation length, X Q , 

is useful in the interpretation of this energy-loss phenomena and is defined by 
Fenyves and Haiman [Ref. 4: p . 25] as: 




Radiation length can be visualized as the distance in which the electron 
energy is reduced to l/e of its initial value and is independent of electron 



energy. In practice it is more convenient to express the quantity in units of 
gm/cm^ (Table 1) in which the calculated values form a smooth curve as a 
function of Z. [Ref. 2: pp. 510-512] 



TABLE 1 

VALUES OF RADIATION LENGTH FOR VARIOUS SUBSTANCES 

Substance Z A Radiation lengths Critical energy 

(gm/cm2) (cm) (Mev) 



Carbon 


6 


12 


44.6 


30.0 


102 


Nitrogen 


7 


14 


39.4 




88.7 


Air 


7.37 


H.78 


37.7 


31.0 xIO 3 


84.2 


Water 


7.23 


14.3 


37. 1 


37. 1 


83.8 


Oxygen 


8 


16 


35.3 




77.7 
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III. ISSUES 



Previous experiments dealing with the experimental verification of 
calculations as they apply to high energy particle beams raised several 
issues. Direct comparisons between experiments and calcualtions have been 
difficult due to differing experimental geometries, varying incident energies, 
and multiple computational computer codes. This experiment will address 
some of these issues. 

A. EXPERIMENTAL GEOMETRY 

The recent experiments by P.F. Cromar [Ref. 1] were conducted in an 
insulated container filled with liquid nitrogen. The use of liquid nitrogen as 
the test medium precluded the immersion of the dosimeters in the 
environmental matrix. The dosimeters were placed in air on the exterior of a 
one inch thick hard foam container. Ideally, direct comparison with 
calculation should be done with the dosimeters placed inside the test medium 
so that the side and back scatter contributions are the same for experiment 
and calculations. 

Separate tanks were required for measurements at different distances 
from the beam entrance. As a result, several independent runs were 
necessary to compile the data at all distances. If the dosimeters could be 
placed directly in the test medium all measurements could be done 
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simultaneously to minimize experimental beam condition variability. For this 
experiment, thermoluminescent dosimeters were sealed in plastic and placed 
directly in water. 

B. ENERGY DEPENDENCE 

The experiment by Cromar was conducted at a single energy of 100 MeV. 
This is roughly the critical energy of liquid nitrogen where there is equal 
contribution from bremsstrauhlung and ionization to the energy loss 
mechanism. The contributions to energy deposition are different from that 
for energy loss. It is possible that the continuous slowing down 
approximation in the computational code may not be adequate at all energies. 

Experiments were conducted at 100 MeV and 20 MeV. The experimental 
results from the 20 MeV run were not returned for inclusion is this paper, but 
results of other published results at varying energies are discussed. 



C. COMPUTATION 

It is not clear from the published data and the calculations of Cromar that 
all versions of the electron transport codes produce the same computational 
results. The apparent disagreement between ETRAN-16 and CYLTRAN needs to 
be investigated. In translating the calculated dose in liquid nitrogen to that 
which might be observed in air, a single geometrical extrapolation was used. 
This may not be adequate as multiple scattering angles and radiation lengths 
in the beam axis and the perpendicular direction interplay. 
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In this experiment, the results of the two-dimensional CYLTRAN code for 
the energy deposition in water have been converted to a one-dimensional 
geometry and compared with the published results for ETRAN. In addition, 
CYLTRAN results for water and liquid nitrogen are scaled and compared to 
the CYLTRAN predictions for air. 



D. DOSIMETRY 

CaF 2 dosimeters have a non-linear response as a function of photon 
energy below 200 keV as shown in Figure 3. 1 . 




Should the radiation field in the vicinity of the electron beam result in a 
large fraction of the deposited energy being due to lower energy photons, 
then the dose measured could be in serious error. The dosimeters have been 
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calibrated with respect to a &0 Co source. The spectrum due to a line electron 
source is not identical to that from 6o Co. 

Dosimeters of CaF 2 and LiF were placed in symmetrical positons for a 

direct comparison of the measured dose. This portion of the experiment 
remains incomplete as the results for the LiF dosimeters were not returned 
due to a faulty set of calibration curves at the reading facility. 
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IV. CALCULATIONS 



A. DESCRIPTION 

The calculations of electron/photon showers were performed using the 
computer program CYLTRAN of the Integrated TIGER Series of Transport 
Codes[Ref. 6], The purpose of the calculation was to determine the off-axis 
dose resulting from discharging a monodirectional beam of electrons into 
several different media. 

The cylindrical geometry of the computer code CYLTRAN is well suited for 
the desired output. A 0.8 cm diameter electron beam of 20,000 particles 
was incident normal to targets of water, air, and liquid nitrogen. Each 
target, as indicated In Figure 4.1, was divided into zones with energy 
deposition being recorded at distances that corresponded to 0.0, 0.25, 0.50, 
1.0, 1.5, and 2.0 radiation lengths in each medium. At each distance, ten 
concentric sub-zones, Figure 4.2, were created to measure the off-axis 
energy deposition. Zone thickness was kept thin to provide semi-discrete 
data points at specific distances of penetration into the medium. However, 
sufficient thickness was necessary to allow for the finite track segments used 
in the computer calculations and to provide for statistical validity of the 
results. 

The output of the code included the target material, zone or sub-zone 
location and mass, energy deposited by primary and knock-on electrons, 
energy deposited by photons, and the total energy deposited within each 
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Figure 4. 1 CYLTRAN Geometry 
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SUB-ZONES 1-10 



Figure 4.2 CYLTRAN Sub-Zones for Water 
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sub-zone. The simulation was run in ten batches of 2000 particles each with 
statistical values being computed after each batch. Each value of energy 
deposition is followed by a two digit number that represents the one-sigma 
uncertainty that results from the Monte-Carlo routine, with low numbers 
indicating narrow bounds of uncertainty. 

At distances well within the targets the statistical uncertainties grew due 
to the fewer number of particles with sufficient energy to penetrate to these 
zones. Decreasing the statistical uncertainties can be accomplished by 
increasing the number of incident particles, however the relationship 
between numbers of particles and uncertainties is not linear. Each 
simulation run of 20,000 particles used about twenty minutes of Cental 
Processor Unit time on the Cray computers at Los Alamos. Increasing the 
number of particles to 45,000 would have reduced the uncertainty by only 
about ten percent but would have increased the. CPU time by a factor of two. 

Since the primary purpose of the computer simulations was to provide 
immediate comparisons with the experimental results, no attempt was made 
to fully optimize the numerical calculations. When reasonable, all default 
values within the codes were used. This led to track segments that were 
larger than ideal at the higher energies near the entrance, while the track 
segments were too small for the less energetic particles well within the 
target. It is believed that for the comparisons desired, the integration of 
the large number of incident particles is adequate to provide meaningful 
results inspite of the less than optimum analysis. 
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Appendix A contains a description of CYLTRAN and sample input data 
streams for the generation of the cross-sections (Table A1) and the 
execution of the program (Table A2). The CYLTRAN output of the energy 
deposition normalized to one incident particle for 100 MeV electrons in water, 
air, liquid nitrogen, and 20 MeV electrons in water is contained in Tables 
A3-A25. 

B. RESULTS 

The total energy deposited in each sub-zone was converted to 
Rads/Coulomb for comparison with the experimental results. The normalized 
dose due to 100 MeV electrons for each sub-zone in water is contained in 
Table B1 of Appendix B. The resultant graphs of the off-axis dose for each 
penetration depth are given in Figures B1 through B7. Appendix’s C, D, and E 
contain similar results for the 100 MeV electrons in air and liquid nitrogen, 
and the 20 MeV electrons in water. 

The curves drawn on the graphs are not intended to necessarily be the 
mathematically best fit approximation for the data although in most cases the 
difference is minimal. Instead, the curves were generated using the 
exponential and polynomial equations of Appendix F, in order to construct the 
three dimensional half-axis surface plots of Figures 4.3 through 4.6 and the 
contour, or isodose, curves of Figures 4.7 through 4.10. The vertical 
error bars result from the indicated one-sigma errors generated by 
CYLTRAN. The plots have been smoothed only in the off-axis direction. 
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Smoothing parallel to the beam axis would have involved a good deal of 
numerical experimentation and iteration and would not have added any 
meaingful interpretation to the results. In spite of the modest smoothing, 
the relatively small number of data points resulted in the rough surface plots 
of Figures 4.7 through 4. 10 and it is felt that no additional smoothing could be 
justified based on these few points. 
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Figure *4.3 : Surface Plot of Dose Due to 100 MeV Electrons in Water 
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Figure 4. 4 : Surface Plot of Dose Due to I 00 lie V Electrons in Air 
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Figure 4.5 : Surface Plot of Dose Due to 100 fleV Electrons in LN 2 
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Figure 4.6 : Surface Plot of Dose Due to 20 MeV Electrons in Water 
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Figure 4.7 : Contour Plot of Dose Due to 100 NeV Electrons in Water 
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Figure *1.8 : Contour Plot of Dose Due to 100 [lev Electrons In Air 
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Contour Plot of Dose Due to 100 MeV Electrons in LN 2 
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OFF-AXIS DISTANCE (CM) 

Figure 4.10: Contour Plot of Dose Due to 20 MeV Electrons in Wate 



V. EXPERIMENTS 



A. DISCUSSION 

All experimental measurements were conducted at the Naval 
Postgraduate School Accelerator Lab using the 100 MeV linear electron 
accelerator. A rectangular container (100 x 46 x 38 cm 3 ) of 4 mm 
polyethylene plastic was used to contain the experimental matrix. This 
allowed for a usable test area of 10 cm on either side of the central axis with 
a minimum of 9 cm of water beyond this to provide a uniform scattering 
medium without side effects. The length of the tank allowed measurements 
out to two radiation lengths in water. 

Calcium fluoride thermoluminescent dosimeters (TLDs) were provided by 
Naval Surface Weapons Center.- The TLDs were wrapped in a single layer of I 
mil aluminum and were enclosed in a thin plastic film to prevent moisture 
from contaminating the crystals. The TLDs were mounted to soft wood 
stretchers at intervals indicated in Figure 5. 1 and Figure 5.2. Wood contains 
hydrogen and carbon and is similar to water in atomic properties. Pairs of 
TLDs (one CaF2 and one LiF) were mounted in symmetric positions. The TLDs 
were generally located at one centimeter spacing, but alternating between 
sides of the central axis, i.e. 1,3,5 cm right and 2,4,6 cm left . This 
served to reduce interference between adjacent dosimeters and was intended 
to provide a means of interpreting the results in case of a slightly off-axis 
beam or drifting of the stretcher from centerline. TLD positions were known 



31 



TLD Positions 



(h 2 o) 



74.0 cm 



_2 ;_0_ R a d_ L_e n^t_h_ 



55.5 cm 



_5_ f? a d_ L_e n^t_h _ 



37 0 cm 



1 .0 R8d Lengt h_ 



1 8.5 cm 



0^ 5 _ Rad _ Length _ 



9 2 cm 



0 25 Rsd Length 



0 0 cm .. 



. P_ _Q . P. 8 J . Len gt_h _ 



tn 

x 

< 



< 

UJ 

a 



A A 

□ □ 



□ 

A 



□ 

A 



A 

□ 



A 

□ 



A 

□ 



□ 

A 



□ 

A 



□ 

A 



AAA 

□ □ □ 



A A A A A 

□ □ □ □ □ 



A A A A A 

□ □ □ □ □ 



A A 

□ □ 



10 8 6 4 2 



□ □ □ 

A A A 



□ □ □ □ □ 

A A A A A 

□ □ □ □ 

A A A A 



□ □ 

A A 



2 4 6 8 10 

(cm) 



A CaF 2 

□ Li F 



Figure 5. 1 TLD Positioning Within Test Tank for 100 Mev Run 
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Figure 5.2 Test Tank Dimensions and TLD Positioning 
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\ 

to within 0.5 cm on the stretcher while the stretcher position was known to 
within one centimeter. 

Pulsed beams of 100 MeV and 20 MeV electrons were normally incident on 
the target at a repetition rate of 60 Hz. The beam passed through a vacuum 
test chamber which was sealed with a five mil aluminum window and then 
through 25 cm of air prior to entering the test tank. 

The total charge delivered was monitored by a Secondary Emission 
Monitor (SEM). A vibrating reed electrometer coupled with a one 
micro-farad capacitor was used to determine the charge delivered to the 
SEM. With an SEM efficiency of .028, the total charge delivered to the matrix 
could be determined. Due to uncertainties in this efficiency factor, it is 
estimated the total charge delivered can only be determined to within five 
percent. Electron energy was controlled by the use of a magnetic deflection 
system and is considered to be sufficiently stable that any fluctuations would 
not be noticable within the limits of this experiment. 

Phosphor screens were placed at both ends of the tank on the central axis 
to align the beam. The TLDs were removed and the beam was passed through 
the empty tank. By remotely observing the phosphor screens, an attempt 
was made to make the beam as parallel as possible while within the tank. 
The beam diameter entering the tank averaged 0.5 cm as indicated by the 
phosphor screens, while beam divergence caused the exiting diameter to 
approach 2.0 cm. 
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Once alignment and focussing were completed the phosphor screens were 
removed, TLDs installed, and the tank filled with water. To prevent 
saturation of the TLDs close to the beam axis and yet still provide a 
meaningful dose on the distant TLDs, three Irradiation segments were 
conducted for each run. Based on estimated dose calculations, selected 
TLDs were removed after each segment to prevent overexposure. 
Normalization of the data to Rads per Coulomb delivered will account for the 
different exposure times. 

/ 

B. RESULTS 

The TLDs were mailed to the Naval Surface Weapons Center, White Oak, 
Maryland for reading. The results of the TLD readings were reported in Rads 
and using the recorded charge delivered to each TLD, the readings were 
normalized to Rads per Coulomb. For the 100 MeV experiment a faulty set of 
calibration curves prevented proper reading of the LiF dosimeters and these 
results were not returned. 

The normalized dose from the CaF 2 dosimeters is contained in Table G! of 

Appendix G. Figures 5.3 through 5.8 are plots of the normalized off-axis 
dose at each length within the water. The vertical error bars are drawn 
based on a ten percent error in the normalized dose . The horizontal error 
bars are the result of the uncertainty in the TLD position on the stretcher 
(±0.5 cm), the strecher location with respect to the central axis of the test 
tank(± 1 .0 cm) , and the skewness of the beam with respect the center line of 
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the tank (±2.0 degrees). The curves drawn on the figures are the 
semi-smoothed results from the CYLTRAN calculations and are included for 
comparison. 

The objective of the experiment was to provide quick but meaningful 
results for comparison with the calculations. As a result, the experimental 
apparatus was kept simple and straight forward. As a consequence, the 
cumulative errors at the extreme distances result in horizontal error bars 
that encompass a large portion of the horizontal axis. However, by placing 
sequential TLDs on alternate sides of the central axis, some simple yet 
justified and reasonable corrections can be applied to the data in order to 
present the true picture. 

Figure 5.3 shows how the true off-axis position of the TLDs varies with 
the degree of skewness in the beam. Since alternate alternate TLDs were 
placed on alternate sides of the matrix x-axis, this skewness results in a 
displacement of the true off-axis distance measured from the centerline of 
the beam. 

Table G2 of Appendix G contains a set of corrected positions for the TLDs 
based on varying degrees of skewness for the beam. A misalignment of only 
two degrees from parallel results in an error of the off-axis position of 2.58 
cm at the furthest distance for which data was recorded. Figures 5.9 
through 5.12 present the results if a uniform correction of 1.7 degrees is 
applied to all the data points. At the inner distances the uncertainty in the 
position of the stretcher is greater than the correction of the beam skewness 
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Figure 5.3 : Effect of Beam Skewness on Off-axis Zero Point 
and the effect of applying the skew correction is not as visable. In all 
cases, the data generally alternates between high and low points of the 
calculated curves, in relation to being on the right side or the left side of the 
central axis. 
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Figure 5.3 : Experimental Results of 100 MeV Electrons in Water-0.0 cm 
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Figure 5.4 : Experimental Results of 100 MeV Electrons in Water-9.2 cm 
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Figure 5.5 : Experimental Results of 100 MeV Electrons in Water- 1 8. 5 cm 
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Figure 5.6 : Experimental Results of 100 MeV Electrons in Water-37.0 cm 
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Figure 5.7 : Experimental Results of 100 MeV Electrons in Water-55.5 cm 
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Figure 5.8 : Experimental Results of 100 MeV Electrons in Water-74.0 cm 
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Figure 5.9 : Experimental Results of 100 MeV Electrons in Water- 18. 5 cm (1.7°) 
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Figure 5. 10 : Experimental Results of 100 MeV Electrons in Water-37.0 cm ( 1 . 
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Figure 5.11: Experimental Results of 100 MeV Electrons in Water-55.5 cm ( 1 .7') 
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Figure 5. 12 : Experimental Results of 100 NeV Electrons in Water-74.0 cm ( 1 . 




VI. CONCLUSIONS 



A. COMPARISON OF COMPUTER CALCULATIONS AND EXPERIMENTAL RESULTS 

Within the limits of this experiment, there is good agreement between the 
measured and calculated normalized dose due to 100 MeV electrons in water. 
The experimental results at 0.0 cm and 9.2 cm were significantly higher than 
predicted. The TLDs at 0.0 cm were intended to measure the dose as a result 
of backscattering. However, the entrance to the test tank was not shielded 
and it is probable that these front TLDs were measuring a significant amount 
of radiation scattered from the aluminum window covering the exit of the 
vacuum chamber. 

The results could be greatly improved with a few simple modifications to 
the TLD matrix. Alternating consecutive TLDs from one side of the central 
axis to the other showed that the beam was not directly on that axis. To adjust 
for this discrepancy it was necessary to try several iterative corrections of 
beam angle to bring the data back in line. If two TLDs had been placed at 
symmetrical positions on either side of the central axis, the two readings 
could be compared , providing a means of calibrating the off-axis positions of 
the TLDs with respect to the true beam position. In addition, use of a rigid 
frame to support the TLDs instead of the floating strechers that were used, 
would allow for much more accurate positioning of the entire array. 
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B. EXTENSIONS OF THE PREDICTIONS FROM ONE MEDIUM TO ANOTHER 

The results were recorded at distances which correspond to fractions of 
the radiation length in the given medium. To extend the results from one 
medium to another, it Is necessary to account for the change in size and mass 
of the detector as the radiation length varies [Ref. 1 : p.34]. If a detector 
subtends the same solid angle and is located at an equivalent distance when 
measured in radiation lengths, the same energy should be deposited in the 
detector. 




R, R 2 



= R j Tan 0 

Area = n r j = .0 l.R | T an 8 ) 



r_ = R ? T an 9 

<■> 

Area = n r\. = .0 ( R ,.,T an 9.) 



Ratio of Areas = 




in (R j Tan 8 ) 
n ( R ,,T an 8 ) 

L 




Figure 6.1 Ratio of the areas of two detectors 



As shown if Figure 6. 1 , the new detector has an area and hence, a mass, 
that is larger by a factor of the square of the ratios of the distance to the 
detector. Since dose is related to the mass of the detector, 

1 Rad = 100 ergs / gm 

the dose at the new detector will be decreased by the same ratio. Water has 
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a radiation length of 37.0 cm and air has a radiation length of 307 meters, so 
for this experiment the ratio is: 



Table 6.1 contains values of the dose in air predicted from the dose in 
water compared to the dose in air predicted by CYLTRAN. At the near 
distances, the results vary by about thirty percent, but this is largely due to 
the high one-sigma errors of the data at these distances. At one radiation 
length and above, the statistical errors in the initial data are much smaller 
and the resulting scaled dose is within ten percent of the dose in air predicted 
by CYLTRAN. Figure 6.2 through 6.6 present curves of the dose in water 
scaled to air, the dose in liquid nitrogen scaled to air, and the dose in air as 
predicted by CYLTRAN. 

The critical energy for water, air, and liquid nitrogen are contained in 
Table 1. Since the experiment and the calculations were conducted at 100 
MeV, close to the critical energies of the three media, the concept of 
radiation length is somewhat nebulous. Water with solid angle effects can be 
scaled to air with consistent results because they have similar Z. It is 
unlikely the results in water would scale to lead, for instance. At lower 
energies, a better length parameter is the range in gm/cm 2 . Figure 6.7 and 
Figure 6.8 contain the off-axis dose in water. Figure 6.9 and Figure 6.10 
contain the results in liquid nitrogen. 
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C . COMPAR ISONS W ITH OTHER PAPERS 



In his CYLTRAN predictions for liquid nitrogen, Cromar included a 0.2 cm 
layer of silicon in the geometry for the energy deposition [REF.l]. This 
special geometry required Individual CYLTRAN runs for each distinct length 
used in the experiment, requiring a great deal of redundant calculations by 
the computer. To avoid this, the calulations in this paper were all done in a 
single run for each different medium. The energy deposition was recorded in 
sub-zones of the medium instead of silicon. To ensure the validity of this 
approach, the off-axis dose in liquid nitrogen predicted by this method is 
compared to the experimental results reported by Cromar. Figures 6.11 
through 6. 13 contain the two sets of data and in general show good agreement 
between the predicted and observed results. 

Cromar's results were also compared to those reported by R.A Lindgren 
[Ref. 7] which looked at the energy depositon in air due to 50 MeV electrons 
using the computer code ETRAN. The differences between the normalized 
doses were significant but several different assumptions were made in the 
two studies and the difference in the incident energy makes a direct 
comparison difficult. 

An additional study of electron transport at energies up to 1000 MeV has 
been reported by S.M. Seltzer and M.J. Berger [Ref. 8]. In their study, the 
computer code ETRAN was used. Since the TIGER series of codes all draw from 
the same source for the Monte-Carlo calculations the results shoud be the 
same as those presented in this paper generated from the CYLTRAN code. 
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Figure 6.H shows the CYLTRAN data plotted with the ETRAN data at incident 
energies of 60 MeV and 125 MeV. The ETRAN curves were derived from 
replicate data of Ref. 8 and plotted at the midpoints of the histogram bins. 
ETRAN is a one-dimensional code and for comparison it was necessary to sum 
the ten radial sub-zones used by CYLTRAN to obtain a single value of energy 
deposition at a given penetration depth. The CYLTRAN zones were only 
0.2 g/cm 2 in thickness so an additional factor of five was included to 
normalize the results to a unit depth. Figure 6.14 shows that once the 
geometries are equalized, there is no significant differences between the two 
computer codes . 

D. ENERGY DEPENDENCE OF THE RESPONSE 

An additional experimental run was conducted at 20 MeV incident energy 
however the results were not returned for inclusion in this paper. For future 
work, it is recommended that the Naval Postgraduate School purchase or 
obtain the use of its own TLD reader and a supply of reusable TLDs. This 
experiment alone used 80 TLDs at each energy. This was was considered the 
minimum necessary to obtain quick, but meaningful results. For a detailed 
study, the number of TLDs could easily run into the hundreds. Having to 
depend on a third party to process this number of TLDs is time consuming and 
imposses a huge burden of the processing facility. It would greatly facilitate 
the studies if the TLDs could be processed at the test facility. 
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E. SUMMARY 



The results presented indicate that the CYLTRAN computer computations 
do accurately predict the experimental results. There does not appear to be 
any dlscreapncy between the one-dimensional computer code ETRAN and the 
two-dimensional computer code CYLTRAN. For materials of similar Z, the 
extension of experimental results in one medium to another appears to be 
valid. 
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TABLE 6.1 



PREDICTIONS FOR AIR BASED ON RADS/COULOMB IN WATER 



Fa ctor=( 30700/37 ) A 2 





Air (Pred) 


One-Sigma 


Air (Cyltran) 


One-Sigma Ratio 


Rad/Lth 


Rads/Coulomb 


Error 


Rads/Coulomb 


Error 



0.00 


7. 9245E+04 


2 


9.280 1 E+04 


3 


1.17 




8. 1257E-0I 


56 


3. 1 579E+00 


91 


3.89 




7.5299E-01 


46 


3.0577E-01 


42 


0.41 




4.6982E-0I 


48 


2.2514E+00 


62 


4.79 




9.9008E-01 


30 


1 .0680E+00 


42 


1.08 




8.3157E-01 


28 


3. I993E-0I 


34 


0.38 




I.2156E+00 


27 


3.8732E-0I 


43 


0.32 




5.5566E-01 


20 


4.8336E-0I 


36 


0.87 




4.0331E-01 


24 


2.4460E-0I 


47 


0.61 




2. 1545E-0I 


27 


1 . 8029E-0 1 
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0.84 



0.25 


8.2900E+04 


3 


9. 5737E+04 


3 


1 .15 




4.0395E+03 


3 


3. 1 162E+03 


6 


0.77 




3.2665E+02 


7 


3.0579E+02 


9 
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8.0173E+01 


13 


8.9421E+01 


18 


1 .12 




3.8481E+01 


16 


3.2072E+01 


28 
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5. 5583E+00 
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40 
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3. 4629E+00 
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• 40 
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3.131 3E+00 


27 


3.2345E+00 


55 


1.03 




2.8272E+00 


52 


3. 2645E+00 


61 


1 .15 



54 



TABLE 6.1 (cont’d) 



2.6881E+04 


3 


2 . 9248E+04 


3 


1.09 


1.2841E+04 


2 


1 .3804E+04 


2 


1.07 


4. 1262E+03 
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4.2394E+03 
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1.03 


1. 1693E+03 
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1. 1804E+03 
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1.01 
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8 
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1 1 
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1 1 
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1 . 10 
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17 
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Figure 6.2 : Predicted Dose in Air (0.0 m) 
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OFF-AXIS DISTANCE (RAD. LENGTH) 
Figure 6. 3 : Predicted Dose in Air (76.8 m) 
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Figure 6.4 : Predicted Dose in Air (153.5 m) 
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Figure 6.5 : Predicted Dose in Air (307.0 m) 
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Figure 6.6 : Predicted Dose in Air (460.5 m) 
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Figure 6.7 : CYLTRAN Dose in Water Due to 100 MeV Electrons 
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Figure 6.8 : CYLTRAN Dose in Water Due to 100 MeV Electrons 
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Figure 6.9 : CYLTRAN Dose in Liquid Nitrogen Due to 100 MeV Electrons 
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Figure 6.10: CYLTRAN Dose in Liquid Nitrogen Due to 100 fleV Electrons 
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DISTANCE FROM BEAM ENTRANCE (G/CM 2 ) 

Figure 6. II : Predicted and Experimental Dose In LN 2 ( 2= 1 .2 g/cmO 
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DISTANCE FROM BEAM ENTRANCE (G/CM 2 ) 

Figure 6.12: Predicted and Experimental Dose In LN 2 (z=4.5 g/cm 
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DISTANCE FROM BEAM ENTRANCE (G/CM 2 ) 

Figure 6. 13 : Predicted and Experimental Dose In LN 2 (z=6.9 g/cm 
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DEPTH Z, G/CM 2 

Figure 6. 14 : Comparison of ETRAN and CYLTRAN Predictions 
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APPENDIX A 



DESCRIPTION OF CYLTRAN AND OUTPUT 

The computer code CYLTRAN is a member of the Integrated TIGER Series of 
Coupled Electron/Photon Monte Carlo Transport Codes developed at Sandia 
National Laboratories by J.A. Halbleib and W.H. Vandevender. CYLTRAN 
employs a fully three-dimensional description of particle trajectories within 
an axisymmetric cylindrical material geometry and finds application in 
problems involving electron and photon beam sources. [Ref. 6: p . 6] 

CYLTRAN is based primarily on the ETRAN model which combines 
microscopic photon transport with a macroscopic random walk for electron 
transport. The code describes the generation and transport of the 
electron/photon shower from the incident energy down to a specified cut-off 
energy for photons and electrons. The code can handle up to ten different 
materials, with the only requirement being that the material geometry be 
cylindrical ly symmetric. [Ref. 6: p . 70] 

The user may specify incident particles to be either electrons or photons 
which may be considered mono-energetic or have a specified energy 
spectrum. The source reference direction can be monodirectional , 
isotropic, or have a specified divergence angle. 
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The default output of CYLTRAN consists of: 

a) Energy and number escape fractions (leakage) for electrons, 
unscattered photons and scattered photons. 

b) Charge and energy deposition profiles. 

c) An explicit statement of energy conservation. 

In addition to the default output, a number of optional outputs may be 
selected through the use of the appropriate keywords. These are: 

a) Escape fractions that are differantial in energy for both electrons 
and scattered photons. 

b) Escape fractions that are differential in angle for both electrons and 
scattered photons. 

c) Coupled energy and angular distributions of escaping electrons and 
scattered photons. 

d) Volume-averaged energy distributions of electron and photon scalar 
fluxes for selected regions of the problem geometry. 

e) Psuedo-pulse-height distributions for selected regions of the 
problem geonetry — for example, those regions corresponding to 
active detector elements. [Ref. 6: pp. 19,34] 

Two basic steps are required for solving any given problem: 

a) Generate the cross sections by running the cross section code. 

b) Run the Monte Carlo code. 

Table A1 is a sample input used for the generations of cross-sections of 
water, standard air (80% nitogen, 20% oxygen), liquid nitrogen, and silicon. 
Table A2 is a sample input for the proper execution of CYLTRAN once the 
cross-sections have been generated. CYLTRAN output for runs of 20,000 
incident electron in water, air, and liquid nitr.ogen is include in Tables A3 
thru Table A25. 
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TABLE A I 



SAMPLE INPUT FOR THE GENERATION OF CROSS SECTIONS 

1 MATERIAL H . 11 1 0 .889 DENSITY 1.0 

2 MATERIAL N .80 0 .20 GAS DENSITY .001226 

3 MATERIAL N DENSITY 0.81 

4 MATERIAL SI 

5 TITLE 100 MeV CROSS SECTIONS FOR H20 AIR N SI 

6 ENERGY 100.0 



72 



TABLE A2 



SAMPLE INPUT FOR THE EXECUTION OF CYLTRAN 



1 

2 

3 

4 
6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 



ECHO 1 
TITLE 

... 100 MeV ENERGY DEPOSITION IN WATER 

******** **************5C)URCE ******************* ******** 

ENERGY 100.0 

* DEFAULT SOURCE PHASE SPACE PARAMETERS 
POSITION 0.0 0.0 0 0.0 

RADIUS 0.4 

* MONODIRECTIONAL Z-AXIS 
DIRECTION 0.0 0.0 0.0 

XXXXXXXXXXXXXXXXXXXXXXQpJ | ONS* * * * * * * * * * * * * * * * * * * * * * * * * * 

HISTORIES 20000 
BATCHES 10 
CUTOFFS 0.05 0.001 
DUMP 

**** ********** *******0UTPUT OPT IONS******************** 
ELECTRON-ESCAPE 
NBINE 2 
NBINT 4 

PHOTON-ESCAPE 
NBINE 2 
NBINT 4 

ELECTRON-FLUX 1 10 
NBINE 2 

PHOTON -FLUX I 10 
NBINE 2 

****** *************»G£OMETRY***** ************ ********* 
GEOMETRY 1 1 

NZ NR ECUT PTCZ 



30 


*ZL 


ZR 


Rl 


RO 




MAT 


IZR 


31 


0.0 


0.2 


0.0 


10.0 


1 


1 1 


10 


32 
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9.2 


0.0 


10.0 


1 
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1 1 


10 
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0.0 


10.0 


1 
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0.0 


10.0 


1 
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10.0 
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TABLE A4 

ENERGY DEPOSITION FROM 100 MeV ELECTRONS AT 9.2 CM OF WATER 
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TABLE A6 

ENERGY DEPOSITION FROM 100 MeV ELECTRONS AT 37.0 CM OF WATER 





to 


to 


to 


tO 


rr 




(XI 


r^ 


to 


(XI 




<\J 


C\J 


(\j 


CXI 


c\j 


00 


(XJ 


(XI 


(XI 


to 


_J 


o 

1 


o 

1 


o 

i 


o 

1 


0 

1 


0 

1 


0 

1 


0 

1 


0 

1 


0 

1 


< 


1 

LlJ 


LlJ 

CO 

00 

in 

(XI 


i 

LU 


1 

UJ 


1 

UJ 


1 

UJ 


1 

LU 


i 

LU 


LU 


1 

LU 


1 — 


cd 


o> 


1 — 


to 


in 


oo 


< — 


i — 


00 


o 


in 


oo 


« — 


CD 


to 


o 


OO 


to 




1 — 


i — 


<J > 


00 


(O 


o 


(XI 


oo 


oo 


< — 




o 


oo 


oo 


(O 


to 


CD 


Tf- 


o 


to 




— 


(XI 


(XI 


(XJ 


(\J 


(XI 


, — 


— 


— 





r^CDintoinoin’— (xi 
a> 
o 



u 


LU 
I ■ 1 


to 


CXI 


CXI 


(XI 


CXI 


to 


to 


to 


to 


to 


(O 

Cl 


LLJ 

CD 

I 


0 

1 


o 

1 


0 

1 


o 

1 


0 

1 


0 

1 


0 

1 


O 

1 


0 

1 


0 

1 


t J 


i 


LU 


LU 


LU 


LU 


LU 


LU 


LU 


LU 


LU 


LU 


c 


<< 


(O 


i — 


OO 


i — 






CXI 


OO 


(XI 


(XI 


a> 


z 


(XJ 


i — 


(D 


O 


in 


O 


to 


(D 


r^ 




TO 


n 


O 


TT 


Tf* 


to 


o 


(XI 


00 


O 


(D 


CD 


3 


< 


Tf 


CXJ 


CXJ 




o 


(D 


TT 


O 




(D 


c 


o 


to 


— 


— 




— 


to 


in 


^r 


(XI 




<1> 
























cz 
























o 




CD 


CD 


C^ 


<D 


< — 


<D 


<D 


(D 


<D 


CD 


o 




CD 


CD 


TT 


(D 


(D 


(D 


CD 


(D 


<D 


C^ 


w 
























TO 

4> 




tT 


tT 


tT 


tT 


Tf 


tJ“ 


Tl“ 


in 




Tt" 


M 


CO 


o 

1 


O 

1 


O 

1 


O 

1 


O 

1 


O 

1 


o 

1 


o 

i 


O 

1 


o 

1 


yfX 


2 : 


LU 


LU 


LU 


LU 


LU 


LU 


LU 


LU 


LU 


LU 


cu 


o 


OO 


OO 


to 


(XJ 


(XJ 


to 


O 


(XI 


(XJ 


(XI 


£ 


1 


to 


00 


tJ“ 


i — 


i — 


to 


(D 




CD 


in 


< — 


XC 


to 


to 


in 


ti- 


O 


to 


0O 


in 


to 


. — 


o 


o 


00 




CD 


(XI 


OO 


« — 


CD 


o 


in 


(XI 


c 


O 


























*— 


CXJ 




to 


Tf 


to 


to 


(XJ 


*— 


in 












1 


| 


1 






1 


1 


> 
























<D 
























H 



























to 


in 


in 


(XI 


tt 


tt 


TT 


r- 


in 


oo 


CO 


to 


to 


(XI 


CXI 


CXJ 


CXI 


CXJ 


CXI 


to 


to 


LU 


O 

i 


o 

i 


o 

i 


o 

i 


0 

1 


o 

1 


O 

1 


o 

i 


o 


0 

1 


cc 


l 

LU 


i 

LU 


i 

LU 


i 

LU 


i 

LU 


1 

LU 


1 

LU 


i 

LU 


i 

LU 


1 

LU 


< 


CD 


in 


in 


in 


to 


— 


to 


— 


*— 




z: 


CD 


oo 


CD 


to 


— 


to 


(XI 


in 


OO 


(XI 




to 


CD 


to 




*— 


Ti- 


to 


00 


00 


TT 


cc 


in 


oo 


in 






to 


to 


o 


Tr 


oo 


CL 


to 


CD 














OO 


in 



LlJ 

O -NK)Trin(DXoo^2 

rxj 



77 





0 




OO 


(XJ 


0 


c^ 


to 


(XJ 


(XJ 


to 




to 


10 


to 


to 


to 


to 


to 


to 


to 


to 


_J 


0 

1 


0 

1 


0 

1 


0 


0 


0 

1 


0 

1 


0 

1 


0 

1 


0 


< 


1 

LU 


LU 

(\J 

CD 


1 

LU 


l 

LU 


1 

LU 


1 

LU 


1 

LU 


1 

LU 


1 

LU 


t 

LU 


1 — 


ro 


^r 


00 


O 


tT 


(XJ 


in 


in 


(\ 


0 


O 


(XJ 


(XJ 


CO 


00 


(XI 


in 


(XJ 


00 


h- 


(XJ 


O 

OO 

• 


OD 


(XJ 


(XJ 


CO 




in 


(XJ 


CD 




to 


IO 


ro 




CD 


CD 


<0 


CD 


CD 




(\j 


tT 


IO 


in 


ti- 


Tf 


to 


to 


(XJ 


, — 



< 

£ 



n 

o 

LO 

in 

in 



< 

CO 

2: 

o 

cr 



r^ 

< 



o 

LU 



1 ■ 1 LU 

s » 

2 * 
o 
o 



a> 

o 



to 

o. 



c 

a> 

“O 



cj> 

C 



a> 

c 

o 

o 



T? 

CD 

N 



H 

o 

oc 

u_ 



O 



CO 

o 

CL 

LU 

O 



TO 

£ 

t_ 

o 

a 



<i> 

51 







in 


OO 


00 


to 


to 


0 


<— 




Ti- 
























O 


ro 


to 


to 


ro 


ro 


ro 


to 


ro 


ro 


ro 


LLJ 

CO 


0 

1 


0 

1 


0 


0 

1 


0 

1 


0 

1 


0 

1 


0 

1 


0 

1 


0 

t 


1 


1 

LU 


1 

LU 


I 

LU 


i 

LU 


1 

LU 


i 

LU 


1 

LU 


i 

LU 


1 

LU 


1 

LU 


< 


O 


«— 


to 


to 


CO 


(XJ 


to 


O 


CD 


TT 


r 


to 


(O 


to 


to 


(XJ 


CD 


(XJ 


r^ 


(\ 


O 


n 


(XI 


Tl" 


(XJ 


(XJ 


in 


(\ 


to 


CD 


CO 


in 


< 


to 


0 


in 


to 


to 


O 


in 


to 


OO 


CD 


0 


(XI 


in 


to 


in 


TT 


in 


TT 


to 


(XJ 


— 




CJ> 


(XJ 


in 


CD 


CD 


CD 


0 


CD 


CD 


0 




CD 


0 


00 


CD 


CD 


CD 


in 


CD 


CD 


r^ 




CO 


ti- 


ti- 


ti- 


tT 


tT 


^r 


in 


in 


in 


CO 


0 

1 


0 


0 

1 


0 

1 


O 

1 


O 

1 


0 

1 


0 

1 


0 

1 


0 

1 


z 


l 

LU 


LU 


1 

LU 


i 

LU 


1 

LU 


J 

LU 


s 

LU 


l 

LU 


1 

LU 


1 

LU 


0 


in 


0 O 


» — 


in 


to 


CD 


OO 




TT 


to 


1 


(XI 


CD 


O 


Ti- 


in 


r^ 


to 


CD 


OO 


TT 




c^ 


IO 


to 


ro 


to 


0 


0 




in 


(XI 


0 


<0 


to 


in 


0 


(XJ 


, — 


CO 


*— 


tT 


«— 


0 
























(XJ 


(XJ 


«— 


to 


(XJ 


* — 


in 


tT 


in 




xc 


t 


t 




l 


l 


l 


l 


1 







CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 



>- 

o 

oc 

LU 

z 

LU 



CO OOOO 

LU OOOO 

— + + + + 

cc LU UJ LU UJ 

< OOOO 

5 ” OOOO 

— OOOO 

cr 0000 

CL . . . . 

OOOO 



OOOOO 
00000 
+ + + + + 
LU LU LU LU LU 
OOOOO 
OOOOO 
OOOOO 
OOOOO 



OOOOO 



LU 

2: O 

O — cum^rinior^oocr)^ 
rxj 



78 



0 . 0000E+00 



TABLE A8 

ENERGY DEPOSITION FROM 100 MeV ELECTRONS AT 74.0 CM OF WATER 
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TABLE A10 

ENERGY DEPOSITION FROM 20 MeV ELECTRONS AT 9.2 CM OF WATER 
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TABLE A 1 2 

ENERGY DEPOSITION FROM 20 MeV ELECTRONS AT 37.0 CM OF WATER 
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ENERGY DEPOSITION FROM 100 MeV ELECTRONS AT 307 M OF AIR 
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ENERGY DEPOSITION FROM 100 MeV ELECTRONS AT 0.0 CM OF LIQUID NITROGEN 
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TABLE A22 

ENERGY DEPOSITION FROM 100 MeV ELECTRONS AT 23.8 CM OF LIQUID NITROGEN 
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TABLE A24 

ENERGY DEPOSITION FROM 100 MeV ELECTRONS AT 71 .5 CM OF LIQUID NITROGEN 
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Figure B7: Unsmoothed Surface Plot of Dose Due to 100 MeV Electrons In Water 
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z axis- 153.5 m 
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Figure D7: Unsmoothed Surface Plot of Dose Due to 100 MeV Electrons In LN 2 
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6. 2832E-01 6.2832E-01 3.4066E-04 5.4218E+07 
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Figure E6: Unsmoothed Surface Plot of Dose Due to 20 MeV Electrons in Water 



APPENDIX F 



EQUATIONS USED FOR CURVE FITTING 
TABLE FI : 100 MEV ELECTRONS IN WATER 



0.0 cm 



0. 5i x i 1.5 


y=1 .71 


x 10 13 


e -1 1 . 49x 


1 . 5< x i 1 0 


y=7. 5 1 6 


x 10 3 


e -0. 106 lx 



9.2 cm 



0. 5i x i 1.5 


y=2 . 59 


x 10 1 ] 


e ~3.022x 


l.5< Xi 2.5 


y=1 .21 1 


x 10 1 1 


e ~2. 51 5x 


2.5< xi 3.5 


y=7.514 


x 10 9 


e -l .405x 


3. 5< Xi 6.5 


y= 1.359 


x 10 9 


e -0.9046x 


6. 5< Xi 7.5 


y=8.291 


x 10 7 


e -0.4732x 


7. 5< Xi 10 


y=5. 1014 


x 10 6 


e -0. 1014x 



18.5 cm 



0.5i x i 1.5 


y=2 . 6675 


x I0 10 


e -0. 7388x 


1 . 5< Xi 2.5 


y=4.8532 


x 10 10 


e -1 . I352x 


2 . 5< Xi 3.. 5 


y=6.6444 


x 1 0 1 0 


e -l . 1 352 x 


3. 5< Xi 4.5 


y=4 . 4876 


x 10>° 


e -1 . 1488x 


4. 5< x i 5.5 


y=6. 7422 


x 10 10 


e -| .2393x 


5. 5< Xi 10 


y=2 . 2332 


x I0 9 


e -0.629x 



37.0 cm 



0. 5i x i 10 



55.5 cm 



74.0 



0. 5i x i10 



0. 5i x i 10 



y= 1 . 848-4 .871 x 1 0 8 x +4 . 807 x 1 0 7 x 2 
1 .748 x 10 6 x 3 



y=4.26 x I0 8 e -0.3299x 



y=2 . 074 x 10° e 



8 o -0.282x 
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TABLE F2 : 100 MEV ELECTRONS IN AIR 



0.0 m 



x i 


12 


y=1 . 587 x 10 7 


e -l.24x 


1 2< x 


i 80 


y=5.972 x 


e -0.0446x 



76.8 m 



4i xi 12 


y=5. 297 


x 10 3 


e -0.4 126x 


I2< x i 20 


y= 1 . 0 1 7 


x 10 5 


e -0.2792x 


20< x i 28 


y=6. 614 


x 10 3 


e -0. 148 lx 


28< x i 36 


y=3.236 


x 10 3 


@-0. 1253x 


36< x i 54 


y=5.50 


x 10 2 


e -0.07838x 


54< x i 80 


y=2 15.3-8. 


. 25x+0. 1066x z -. 000458: 



153.5 m 

4i x i 12 
1 2< x i 20 
20< x i 28 
28< x i 36 
36< x i 44 
44< x i 80 



y=4.262 x 10 4 
y=8. 1 107 x 10 4 
y- 1.0363 x 10 5 
y=8. 5209 x 10 4 
y=5. 6628 x 10 4 
y=4. 1004 x 10 3 



e -0.0905x 
e -0. 1422x 
e ~0. 1 540x 
e -0. I473x 
e -0. 1364x 
e -0.0789x 



307.0 m 



4i x i 80 



y=2 . 802 x 10 3 



e-0. 04468x 



460.5 m 



4i x i80 



y=6. 90 x 10 2 e -0.04427x 
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TABLE F3 : 100 MeV Electrons in Liquid Nitrogen 



0.0 cm 

0.5i x i 1.5 
1 . 5< x i 10 



11.9 cm 

0.5i x i 2.5 
2. 5< x i 3. 5 
3. 5< x i 4. 5 
4. 5< xi 10 



23.8 cm 

0.5i x i 1.5 
1 .5< x i 2.5 
2. 5< x i 6. 5 
6. 5< x i 10 



47.7 cm 

0.5i x i 10 



71.5 cm 

0. 5i x 1 10 



95.4 

0.5i x i 10 



y=6. 1 66 


x 10*2 


e - 1 1 . 49x 


y=4. 198 


x 10 6 


e -0. 106 lx 



y=1 .593 


x 10 1 1 


e -3.022x 


y=2. 162 


x 10 10 


e -2.5l5x 


y=5. 198 


x 10 9 


e -l .405x 


y=5. 525 


x 10 8 


e -0.9046x 



y=1 .369 
y=2.397 
y=3.873 
y- 1.853 


x 10 10 
x 10 10 
x 1 0 1 0 
x 10 10 


e ~0.7388x 
e ~1 . 1352x 
e -1 . 1352x 
e -1 . 1488X 


y= 1 . 17 


x 10 9 


e ~0.258x 


y=2.61 


x 10 8 


e -0.2817x 


y=1 .221 


x 10 8 


e ~0.23l9x 
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TABLE F4 : 20 MeV Electrons in Water 



0.0 cm 



9.2 cm 



18.5 cm 



37.0 cm 



55.5 cm 



0. 5i x i 1.5 
1 . 5< x i 1 0 



y=5.4289 x 10 12 e _9 - 198,x 

y= 1.91 99 x 10 7 e "0.4918x 



0.5i x i 1.5 


y=4. 1483 


x 10 9 


e -0.3287x 


1 . 5< x i 2 . 5 


y=8.3790 


x 10 9 


e -0.7938x 


2. 5< x i 3.5 


y=2. 1945 


x 10 10 


e -1 . 1825x 


3. 5< x i 4.5 


y=1 .7157 


x 10 n 


e ~1 . 77x 


4. 5< x i 6. 5 


y=1 .6720 


x 10 10 


e ~1 . 262x 


6. 5< Xi 7.5 


y=6 . 9 1 53 


x 10 8 


e ~0. 7654x 


7. 5< xi 8.5 


y=5. 8735 


x 10 8 


e ~0.7436x 


8. 5< x i 10 


y=5.0129 


x 10 6 


e -0. 1 832x 


0.5i x i 10 


y=6. 5584 


x 10 7 - 2. 


1939 x 10 7 x + 



2.6384 x I0 6 x 2 - 1.0869 x 10 5 x : 



0.5i x i 10 


y=7. 572 


x 10 6 


g -0. 209x 


0. 5i. x i.10 


y=8. 286 


x 10 6 


e -0.3157> 
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APPENDIX 6 



EXPERIMENTAL RESULTS 

Table Gl contains the results of the 100 MeV experimental run. The 
uncertainties involved include the TLD position on the stretcher (± 0.5 cm), 
the strecher location with respect to the central axis of the test tank (± 1 .0 
cm), and the skewness of the beam with respect the center line of the tank 
(±2.0 degrees). Uncertainty in the charge delivered is five percent and 
uncertainty in the normalized dose is ten percent. 
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TLD 

Number 

I 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

13 



TABLE G1 



Experimental Run # 1: 100 Mev Electrons in Water 



Test day - 28 JUNE 1985 

Test Medium - Water 

Beam Diameter -0.5 cm 

(0.4 cm above the horizontal axis) 

SEM Efficiency - 0.028 

Capacitor - IxlO -6 Farad 

TLD Type - CaF 2 



z- axis 


x-axis 


Rads 


Coulomb 


(cm) 


(cm) 







0.0 


2.0 


35.6702 


1 .03E-06 


0.0 


4.0 


15.8316 


1 .03E-06 


0.0 


6.0 


64.0929 


1 .36E-05 


0.0 


8.0 


49.6843 


1 .36E-05 



9.2 


1 .0 


2348.76 


1 .03E-06 


9.2 


2.0 


631.758 


1.03E-06 


9.2 


3.0 


183. 1 13 


1 .03E-06 


9.2 


4.0 


78.6993 


1 . 03E-06 


9.2 


5.0 


42.9735 


1 .03E-06 


9.2 


6.0 


22.6433 


1 .03E-06 


9.2 


7.0 


18.3424 


1 .03E-06 


9.2 


8.0 


111.069 


1 .36E-05 


9.2 


9.0 


111.267 


1 .36E-05 



Rads/Coul 



3.46E+07 
1 .54E+07 
4.71E+06 
3.65E+06 



2.28E+09 
6. 13E+08 
1 .78E+08 
7.64E+07 
4. 17E+07 
2.20E+07 
1 .78E+07 
8. 17E+06 
8. 18E+06 
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TABLE 61 (cont'd) 



TLD 


z- axis x-axis 


Rads 


Coulomb 


Rads/Coul 


Number 


(cm) (cm) 









14 


18.5 


1.0 


12188. 1 


1 .03E-06 


1 . 18E+I0 


15 


18.5 


2.0 


4085.73 


1.03E-06 


3.97E+09 


16 


18.5 


3.0 


2407.99 


1.03E-06 


2.34E+09 


17 


18.5 


4.0 


516.-688 


1.03E-06 


5.02E+08 


18 


18.5 


5.0 


297.201 


1 .03E-06 


2.89E+08 


19 


18.5 


6.0 


97.8459 


1 . 03E-06 


9.50E+07 


20 


18.5 


7.0 


453.528 


1 .36E-05 


3.33E+07 


21 


18.5 


8.0 


340.428 


1 .36E-05 


2.50E+07 


22 


18.5 


9.0 


194.366 


1 .36E-05 


1.43E+07 


23 


18.5 


10.0 


156.074 


1.36E-05 


1. 15E+07 



24 


37.0 


1 .0 


1282.73 


1.03E-06 


1 .25E+09 


25 


37.0 


2.0 


1783.29 


1 .03E-06 


1.73E+09 


26 


37.0 


4.0 


3612.04 


1 .36E-05 


2.66E+08 


27 


37.0 


6.0 


5309.53 


1 .36E-05 


3.90E+08 


28 


37.0 


8.0 


4737.13 


1 .02E-04 


4.64E+07 


29 


37.0 


10.0 


7204.42 


1 .02E-04 


7.06E+07 



30 


55.5 


1 .0 


5743.79 


1 .36E-05 


4.22E+08 


31 


55.5 


2.0 


1 189.37 


1 .36E-05 


8.75E+07 


32 


55.5 


4.0 


13359.7 


1 .02E-04 


1 .31E+08 


33 


55.5 


6.0 


3039.64 


1 .02E-04 


2.98E+07 


34 


55.5 


8.0 


4717.41 


1 .02E-04 


4.62E+07 


35 


55.5 


10.0 


1047.85 


1 .02E-04 


1 .03E+07 



36 


74.0 


2.0 


2625.15 


1.36E-05 


1 .93E+08 


37 


74.0 


4.0 


2625.15 


1 .02E-04 


2.57E+07 


38 


74.0 


6.0 


7757. 14 


1 .02E-04 


7.61E+07 


39 


74.0 


8.0 


1181.88 


1 .02E-04 


1 . 16E+07 



141 



INITIAL DISTRIBUTION LIST 



No. Copies 

1. Defence Technical Information Center 2 

Cameron Station 

Alexandria, Virginia 22304-6145 

2. Library, Code 0142 2 

Naval Postgraduate School 

Monterey, California 93943-5100 

3. Professor F.R. Buskirk, Code 61 Bs 5 

Department of Physics 

Naval Postgraduate School 
Monterey, California 93943-5100 

4. Xavier Maruyama, Rm. B108, Bldg 245 5 

Center fo Radiation Research 

National Bureau of Standards 
Gaithersburg, Maryland 20899 

5. Professor J.R. Neighbors, Code 61 Nb 2 

Department of Physics 

Naval Postgraduate School 
Monterey, California 93943-5100 



6. Mr. Donald Snyder, Code 61 Ds 1 

Department of Physics 

Naval Postgraduate School 
Monterey, California 93943-5100 

7. Joe Mack 1 

Code M-4 

Los Alamos National Laboratory 
Los Alamos, New Mexico 87545 

8. Gene Nolting 2 

Naval Surface Weapons Center 

Code H23 

White Oak, Maryland 20810 



142 



9. Ms. Louise Miles I 

Naval Surface Weapons Center 

Code H23 

White Oak, Maryland 20810 

10. Or. Ken Struve 1 

P.O. Box 808 

Lawrence Livermore National Laboratory 
Livermore, California 94550 

11. Dr. R. A. Lindgren I 

Department of Physics 

University of Massachussets 
Amherst, Massachussets 01033 

12. Lt R.D. Fitzpatrick 2 

647 Tewkesbury Ln 

Severna Park, Maryland 21 146 



143 



) 



/ / ) 



Ti iesi s 

F4837 




Study 0 f f 
rad iati on exn f Xls 
from relat- • X ?° Sure 

thr »^h matter. erS ^ ng 



Thesis 

F4837 Fitzpatrick 
c.l Study of off-axis 

radiation exposure 
from relativistic 
electrons traversing 
through matter. 




thesF4837 

Study of off-axis radiation exposure fro 



3 2768 000 64959 4 

DUDLEY KNOX LIBRARY 








